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^ I Abstract. In this article we investigate the principle and properties of a vertical 

passive seismic noise attenuator conceived for ground based gravitational wave 
interferometers. This mechanical attenuator based on a particular geometry of 
cantilever blades called monolithic geometric anti springs (MGAS) permits the design 

CsJ i of mechanical harmonic oscillators with very low resonant frequency (below lOOmHz). 

^ \ Here we address the theoretical description of the mechanical device, focusing on 

the most important quantities for the low frequency regime, on the distribution of 

f^ . internal stresses, and on the thermal stability. 

^«0 I In order to obtain physical insight of the attenuator peculiarities, we devise 

some simplified models, rather than use the brute force of finite element analysis. 
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^^ , Those models have been used to optimize the design of a seismic attenuation system 

"^ ' prototype for LIGO advanced configurations and for the next generation of the TAMA 

^^. interferometer. 
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o3 ■ 1. Introduction 



A fundamental issue in designing interferometric gravitational wave detectors is the 
isolation of the test masses (heavy mirrors) from environmental noise sources, such 
as seismic noise. It is well known that a very effective way to reduce this type of 
perturbation is to use mechanical harmonic oscillators, which possess the important 
feature of attenuating the seismic noise above their resonant frequency. 

One of the major difficulties in designing this kind of system is to achieve vertical 
attenuation performance comparable to that obtainable in the horizontal direction. In 
long baseline interferometers, this is a significant issue because the vertical motion 
is coupled to the direction of sensitivity (horizontal direction) through mechanical 
imperfections, optic wedges, the earth curvature, and mirrors' orientation. 
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Figure 1. Schematic view of a complete MGAS filter. 

A possible solution to the problem, which has been implemented in the so-called 
super-attenuators of the VIRGO detector [Caron et al 1997], is based on the magnetic 
anti-spring concept [Beccaria et al 1997]. 

A considerable improvement to the magnetic anti-spring is the geometric anti-spring 
concept (GAS) [Bertolini et al 1999, Cella et al 2002], which allows significant increases 
in the vertical attenuation and thermal stability of a seismic attenuation system (seismic 
filter). 

The monolithic geometric anti-spring (MGAS) concept introduced by one of the 
authors is essentially an evolution of the GAS concept capable to significantly improve 
the GAS system performance. 

A seismic attenuation system based on the MGAS concept is better than the 
equivalent GAS system mainly because of its simplicity. In fact, this characteristic allows 
removal of several mechanical parts that introduce unwanted low frequency resonances 
responsible for reducing the filtering efficiency on all the degrees of freedom. 

In the next sections we address the theoretical description of the MGAS, and some 
of the main issues related to the design of a seismic filter. 

2. Low Frequency One Dimensional Model 



An MGAS configuration is made of several quasi-triangular blades radially disposed and 
connected together at their vertices (see Figure HJ. 

Taking advantage of this radial symmetry, the MGAS physics can be studied 
considering just a single blade. Moreover, because of the blade's symmetry along a 
vertical plane (called Oxy as shown Figure EJ, and the typical blade dimensions (thin 
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blade approximation), we can reduce the dynamics to one of a simple single degree of 
freedom system. 
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Figure 2. Reference Frame for the MGAS study. The curvilinear coordinate I is 
measured from the point O and the component F^, and Fy represent the external 
forces appHed to the blade 



The blade's base is clamped to a massive structure and subjected to an external 
force on the tip due to the load and to the other radially arranged blades. The horizontal 
component compresses and bends the blade longitudinally. 

In the low frequency regime, below the internal mode frequencies, the dynamics is 
dominated by the clamping structure mass (much greater than the blade's mass, which 
we can neglect) and by the load. 

The blade has a constant thickness d much smaller than the two other dimensions, 
a length L and a variable transverse width w{l), with / G [0,L]. The variable / is the 
curvilinear coordinate following the intersection of the blade's neutral surface with the 
plane Oxy, i.e. the neutral axis of the blade. It is also convenient to define 9{l) as the 
angle between the tangent of the curvilinear coordinate and the x axis. 

Under these approximations, we can model the thin blade as a massless elastic line 
with a potential energy [Landau, Lifshitz 1959] 
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Fx and F^, are respectively the horizontal and vertical external forces applied to the 
blade's tip, and I is the transverse moment of inertia of the blade, / = w{l)d^ /12. The 
mechanical characteristics of the blade material enter just through the Young's modulus 
E. 

It is convenient to rewrite the potential energy in the following way 
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where 
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and 7(x) = w{0)/w{x) is a normalized blade's shape function, ^ = l/L. The quantities 
Gx-iGy are dimensionless parameters proportional to the horizontal compression and to 
the vertical load: 

12L2 

^' = E¥M^f- ^ = "'^- ^') 

Writing U in such a form we make clear that the nonlinear behavior of the blade 
(the U term) is parameterized by the dimensionless parameters Gi, which completely 
defines the "working point" of the system for a given normalized shape of the blade 7(0- 
Every physical quantity we are interested in will be indeed written as a product of two 
terms. The first is an appropriate dimensional scale factor, a simple combination of the 
blade's parameters. The second is a function of the G, only (for a given shape of the 
blade), and must be numerically evaluated. This means that the solution corresponding 
to a particular pair {Gx, Gy) can be used to describe a large class of different physical 
systems, obtained one from the other by a transformation of the parameters which leaves 
Gx and Gy both unchanged. For example, the vertical position Y of the blade's tip can 
be written as 

Y = L f\m9{0d^ = Ly{G,). (5) 

Scaling the blade length L by a factor A, and multiplying the blade's thickness d hy a 
factor A^/^ , we keep the working point unchanged. This is very useful in order to plan 
the integration of the blade into a new system [Bertolini 2003]. 

2.1. Integration of the Equation of Statics and General Properties 

Computing the variation 6U of the action ^ with respect to 9 to minimize the potential 
energy, we obtain the equation of statics for the physical system in the presence of an 
assigned external load Gj. Then, rewriting 6U as a system of first order differential 
equations we get 

^ = G. sin0(O-G, cos0(O, (6) 

dO ,,, ,_, 

^=7«)P. (7) 

with boundary conditions 

0(0) = 00, 0(1) = 01 (8) 

Because this is a boundary value problem, in general the solution uniqueness cannot 
be guaranteed. Different equilibrium configurations can coexist for the same values of 
boundary conditions and parameters. This can be easily seen in the particular case 
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of constant blade width. In this case, 7 = 1 and the system admits a ^-independent 
quantity A that can be written as 



where we set /3 = O — ijj, G = JGI + G^, r = v^^ and tan?/) = Gy/Gx- This is formally 
equivalent to the energy of a simple pendulum and we will discuss now the behavior of 
our system from a qualitative point of view, looking at the corresponding phase diagram 
depicted in Figure 01 In [Appendix B we give some details about the analytical solution 
of this problem. 




Figure 3. Evolution of the blade's angle in the phase space, for a blade of constant 
width. The phase space is equivalent to the one of the simple pendulum: closed 
trajectories correspond to a bounded variation of the blade's angle, open ones to blade's 
configurations with turns. Different solutions which satisfy the boundary condition for 
/3(0) = /3o and /3(vG) = /3i are depicted: solution Ti starts in Si and ends in Ei. We 
assume periodic boundaries in /?: Fi starts in 5'i and after a complete 27r cycle moves 
to El. See Section WA\ and p^ppendix B| for a full discussion. 



Our formal "pendulum" is allowed to evolve only in a well defined interval of the 
"time" parameter, that is r G [0, vG]- During this interval the system must move from 
a starting point on the [3 = [3q line to a final point with /3 = /5i. We can freely choose 
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our starting point on the line (3 = (5q. Different possibilities correspond of course to 
different values of A. 

Suppose that we start on the boundary of the oscillatory region (the ^o point in 
Figure IS]). Then, we find the first solution Fq when G is large enough to allow the end 
point to be in Eq. This is a very simple configuration, with a monotonically increasing 
value of 6'. If we now keep G constant and increase p at the starting point we obtain 
a motion which is faster and faster, and a longer and longer trajectory. When p is 
big enough, the trajectory will end at /? = /?i + 27r, describing a new solution Fi with 
the same parameter G as the old one. Clearly, with this method we can generate an 
infinite set of solutions F^, which differ in the number of times the blade's angle makes 
a complete turn. Configurations with turns can potentially give good working points 
(see for example [Winterfiood and Blair 1998]), but we are not interested in these in 
this work. 

Another possibility is the existence of "undulated" configurations. Inside the 
oscillatory region of the phase diagram we can construct solutions with the same (large 
enough) value of G, which can be classified according to the number of zero curvature 
points (crossings of the /3 = axis in the phase diagram) . To make this plausible we can 
consider a solution with some starting point S^ (see Figure [3) and a very large G, such 
that it makes N turns in the phase diagram. When we move S^ towards the boundary 
of the oscillatory region, we should find at least 2N solutions. The reason is that on 
the boundary the solution is without crossings, because the time required to approach 
the point /3 = is infinite. So the trajectory must "unwind" when we move 5*3, and its 
end point must cross the (3 = jSi line at least 2N times. Connected with the existence 
of "undulated" configurations is the possibility of buckling-type phenomenaff. As we 
will see, this is in fact the key to understanding the behavior of our system, as will be 
discussed extensively in the following section. 

2.1.1. Effective Potential To study the solutions of Equations (jH)) and (|2l) as a function 
of the parameters, it is convenient to apply the so-called continuation method. We 
use the AUTO2000 code, which implements several algorithms for continuation of 
differential and algebraic problems. Details can be found in the user's manual [Doedel 
et al 2000]. 

The basic strategy consists in the exploration of the solution space when smoothly 
changing the values of the coordinates of the blade's tip, the horizontal position 

x = ^= j\ose{i)di (10) 

L Jo 
or the vertical position defined in Eq. Q. 

Using the continuation method we constructed the plot in Figure HI which is a 
graphic representation of Gy{x,y). In the explored region, it appears that this is a 
well defined single- valued function as suggested by physical intuition. In fact, numerical 

ft See [Winterfiood et al 2002] for a different approach to vertical seismic attenuation, based on a 
buckling behavior. 
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Blade tip horizontal position normalized to blade length X/L 



Figure 4. The normalized vertical position y versus normalized horizontal position 
X of the blade's tip, for different values of the vertical dimensionless load Gy. For a 
generic value Gy there is a folding singularity where a couple of solutions vanish, and a 
stable branch. In the special case Gy = G* (here G* ~ 1.6895) we obtain a two sided 
bifurcation with a turning point. This is depicted schematically in the lower square 
(bold lines), emphasizing its difference from the simpler pitchfork singularity depicted 
in the upper square. See Subsection 12. 1.61 for a full discussion of these diagrams. 



investigation using AUTO2000 gives no evidence of bifurcations. We can repeat a similar 
study to obtain a representation of the function Gx{x,y) in the same region. We do 
not show the relative diagrams because they do not present peculiar features, but the 
final conclusions are that both G^ and Gy are regular and single-valued functions of 
the blade's tip coordinate. We can expand them in small variations of the parameters 
around the chosen working point 
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where S^ = x — Xq and 6y = y ~ y^. As Gj are the components of the external force 
acting on the system, it is natural to ask if they can be written as the partial derivatives 
of an effective potential function 

dUeff{x,y) 



Giix.y) 



dxi 



(12) 



On physical grounds, the answer should be positive. In fact, this potential energy 
represents the energy stored in the system by the external work, and because there is 
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no dissipation in our model, it should be a well defined function of the position. If 
Ueff were multivalued, the difference between two of the energy values could only be 
interpreted as the difference between the energy stored in two different configurations 
with the same x,y values. But in this case a closed path in the (x, y) space crossing a 
bifurcation point should exist. If Ueff exists the integrability conditions tell us that the 
coefficients K... are completely symmetric in their indices. 

A comment about the general validity of these results is in order. As we discussed, 
it is possible that many solutions with the same values of x and y exist, so Equations (fTTll 
should be understood as a perturbative expansion around some well defined blade 
configuration. In practice, this all we need. With this restriction all should work well, 
unless we try to expand around a bifurcation point (that is, a point where the number 
of solutions changes). 

All the quantities we are interested in can be recovered from the coefficients K..., 
which are partial derivatives of the effective potential evaluated at the working point. 
It is quite important to find a good method to monitor these quantities during the 
continuation procedure. This can be done by introducing auxiliary functions, governed 
by appropriate differential equations that must be added to the basic relations ^, ^ 
and ((HI). The method is explained in detail in [Appendix A| here we stress only that 
for a given order in the expansion of the effective potential we obtain a finite system of 
differential equations to solve. 

As a final observation, we recall that the assumption of the existence of an effective 
potential is not mandatory in order to calculate the coefficients K.... In fact, the 
integrability conditions translate into a set of nontrivial relations between the auxiliary 
functions which can be used to test or simplify the numerical code (see [Appendix A| ). 

We will give other details about the continuation method in Subsections 12.1.21 and 
I2.1.3[ explaining its application to the search for configuration with a low vertical 
stiffness. For more information, good starting points are [Doedel et al 1991, Doedel 
et al 1991b] and references therein. 

2.1.2. Vertical Position of the Blade's Tip To give a simple example of application of 
the continuation method, we will study the normalized blade tip's vertical position. 

Here and in the following we will present results for a particular configuration with 
6q = 7r/4, 9l = —it/6 and 7~^(0 = ci + C2Cos/5^ + cssin/?^, with ci = —0.377, 
C2 = 1.377, C3 = 0.195 and P = 1.361. The particular shape of the blade was 
chosen to compare with measurements taken on a prototype which will be published 
elsewhere [Sannibale et al 2004]. In Subsection 12.1.61 we will comment on the general 
validity of our conclusions. 

If we choose as independent parameters the vertical load Gy and the tip's horizontal 
position X, we can construct the diagram in Figure [HI There the possible equilibrium 
values of y are plotted as a function of the vertical load, for different values of x, in an 
"interesting" region. 

It is evident that the curves can be grouped in two different families, one made of 
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monotonically increasing functions of Gy (continuous line curves, corresponding to large 
values of x) and the other made of curves showing a couple of folding singularities 
(dashed line curves, folds marked with circles). The boundary between these two 
different families corresponds to x = x* ~ 0.9027. 

A folding singularity can be detected during the continuation procedure performed 
by AUTO2000. After the detection it is possible to add a second free parameter and 
to continue the fold. This is very convenient to determine the locus of all the folding 
singularities. 

When X > X* there is only a single equilibrium position yo for each value of the 
vertical load, and as dy/dGy > we can conclude also that the equilibrium is stable. 
If we adiabatically reduce x, the blade's y will remain at the assigned vertical load: yo 
will increase if Gy > G* = —1.6895 and it will decrease otherwise. 

When X < X* it becomes possible to have more than one equilibrium configuration: 
in that case two new solutions yi and y2 appear, and by checking the sign of dy/dGy 
it is easy to see that only the largest and the smallest among yo, yi, y2 are stable. The 
system can switch between these if driven by a large enough external perturbation, and 
we conclude that the region x < x* is associated with a bistability. The presence of 
this bistability can be described also by observing that the system can show a kind 
of hysteresis behavior. This is also illustrated in Figure El suppose we prepare the 
system in a configuration which corresponds to the point (a). If we start to increase the 
load adiabatically, we observe that the vertical blade tip's position increases smoothly, 
moving along the curve, until the fold singularity (6) is reached. At that point the system 
cannot respond anymore in a smooth way to a further increase of the load: it will switch 
abruptly to the only available equilibrium position (c). Of course the switch between 
(6) and (c) points is a dynamical process and depends on the system inertia. It cannot 
be described without solving the nonlinear equation of motion including dissipation 
mechanisms. Our static analysis only shows that when the equilibrium is reached the 
system must be settled in (c). 

If now we reduce the load, the system moves smoothly along the new curve until 
it reaches the folding point (d). A further load reduction triggers a switch to the point 
(a), completing the cycle. 

2.1.3. Vertical Resonance Frequency A very important parameter of an MGASF filter 
is the angular frequency oOy of the fundamental vertical mode. For symmetry reasons it 
is evident that during the filter's vertical motion the horizontal position of the blade's 
tip is fixed, so the equations which define our problem are ^ and ^, supplemented 
by the two fixed boundary conditions ^ and by the geometrical constraint (fTO|l . The 
external parameters are Gy and x, while Gx is fixed by the constraint. 
It is easy to see that in the frequency domain, and for low frequency, 

L J_(dFy\ _J_(dGy\ _Ky^ 

g'^y- Fy\dy)^-Gy\dy),^- Gy- ^''^ 
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Figure 5. The normalized vertical position y of the blade's tip versus the 
dimensionless vertical load parameter Gy (see Eq. ^ . Each curve corresponds to a 
different fixed horizontal position x of the blade's tip (the x value is shown above 
each single curve). Two different classes of curves can be distinguished. For a curve 
in the first class (depicted with a continuous line) the vertical position is always an 
increasing function of the load, so that each point defines a stable equihbrium point for 
the system. For a curve in the second class (depicted with a dashed line) there is also a 
region, where y decreases with Gy, of unstable equilibrium points. This unstable region 
is separated from the stable one by a couple of folding singularities. The continous 
folded curve, which corresponds to x = 0.9027 in this example, separates the two 
regions. The oriented closed path shows a typical hysteresis path, and is discussed in 
Subsection l2.1.2l Note that the two stable equilibrium positions in the bistable region 
correspond to different vertical positions of the blade's tip. Boundary conditions are 
always 6*0 = 7r/4 and 9l — — 7r/6. 



^y^^f^y (16) 
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In other words, the MGAS filter is indeed formally equivalent to a simple pendulum 
with eff'ective length 

Leff = LGyj^^-y. (14) 

\ dy Jx 

For the optimization of filter performance, we are interested in obtaining the 

smallest possible value of Uy. As already stated, this can be obtained by simply 

increasing Lejf while keeping the working point Gi fixed. A possible transformation 

that performs this task is a rescaling of the blade length by a factor A 

F IF 

T ^ XT ""^y ) ""'^ (^^] 

' EdM^) A2 EdM^) ^ ' 

which changes the frequency accordingly with 

1 

7a' 

The vertical frequency of the blade decreases only with the square root of the blade 
length. A long blade is quite difficult to integrate into a real system, and its internal 
resonances start to reduce attenuation performance at low frequency. For these reasons 
we are interested in working points where Lg// diverges to infinity. These are the points 
of the parameters space where the linearized response of the system is large, or in other 
words where the coefficient Kyy of the effective potential goes to zero. 

To describe the typical situation we are looking for, let's re-analyze FigureEl From 
Equation (fT^ it is evident that the effective length of the system is inversely proportional 
to the slope of a curve at the considered equilibrium point. This means that in a folding 
singularity, Gy = G^J°^'^\x), the effective length diverges and so Uy -^ 0. Another point 
where Lg// — i> oo is on the boundary curve x = x* , for the particular value Gy = G* 
which corresponds to a vertical slope. 

Setting the parameters of our system "near" one of these critical points, we can 
obtain a vertical resonance frequency as small as we want in principle. However, there is 
a qualitative difference between the critical points on a folding singularity (x, G^^°'''^\x)) 
and the critical point on the boundary curve {x*,Gy). 

This difference can be fully appreciated by looking at Figure IHl Here the square 
of the vertical resonance frequency in g/L units is plotted versus the vertical load 
parameter, for different x values. In this picture it is easy to distinguish between the 
regions of stable equilibrium points (Uy > 0) and the unstable one {Uy < 0). An 
intersection between a curve and the axis ujy = corresponds to a folding singularity in 
Figure El 

We can choose the working point near this intersection, in the stable region: however 
in this case a load fiuctuation can easily drive the system into the unstable region, and 
this is dangerous because we are in a bistable regime. On the contrary the critical point 
on the boundary curve corresponds in the plot to the cusp on the Uy = axis. If we 
choose our working point near this, on a fully stable curve x = x* +6x, Sx > 0, it cannot 
be driven by a load fiuctuation into the unstable region. Of course, a fiuctuation of x 
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Nonnalized vertical load G (dimensionless) 



Figure 6. Vertical frequency dependence versus load, near the minimum frequency 
working point. The vertical axis is the squared angular frequency of the vertical 
oscillation, in the low frequency approximation, expressed in units g/L. Each curve 
is labeled by its x value. Dashed curves correspond to a; > x* , continuous ones to 
X < X* . The bold curve, with a cusp at ujy = 0, is the boundary between the folded 
and unfolded region. The same hysteresis cycle of Figure is plotted, with the same 
labels {a, b, c, d). In the upper left corner we plot a detailed view of the region around 
the frequency minimum. 
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(which can be generated by a horizontal motion of the blade's tip) can drive our system 
into the unstable region, but in principle this is less dangerous because we don't have 
bistability in this region and the system should fall back to the original point after the 
fluctuation. This does not means that there are no problems. The main issue here is the 
effect of the nonlinearities that unavoidably appear in systems with stiffness approaching 
to zero. This is discussed in section 12.1.41 

We can parameterize Figure IHl with the K*, coefficients. Looking at Figure El we 
see that in the region of interest Gy is accurately described by a cubic function of y 
with coefficients that depend on x. If we take as reference point (x*,^*) and neglect 
corrections in 6x we can write 

1 3 

Gy = -Kyyyy6y + G^ ( 1 7) 

and for the frequency 

-^y - — T^ (18) 

9 f^y 

eliminating 6y we finally get near the working point 

2/3 

(19) 



L K^/^ k{Gy-G 

^ , ,2 _ ^^yyyy L V ^ 

— '^ — ;:: 
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where the cusp singularity is evident. Using the same procedure we can obtain a 
general relation, more involved, which reproduces quite accurately Figure El for small 
enough 5^ 7^ 0. 

It is interesting to add some considerations about the bistable regime, drawing 
the same hysteresis cycle we discussed in Subsection 12.1.21 Looking at Figure El it is 
evident that jumps are located at the edge of the unstable region, and the two stable 
equilibrium points correspond not only to a different blade tip vertical position, but also 
to a different value of the vertical resonance frequency. An exception is, for example, 
the point D. 

2.1.4. Nonlinearities In Section f2.1.31 we saw that from the point of view of system 
stability the critical point on the boundary curve {x*,G*) should be preferred over 
the critical point near a folding singularity. There is another good reason of choosing 
{x*, Gy). The fundamental peculiarity of the MGAS is that the variation of the stiffness 
with the position is in principle not negligible. This means that around the working 
point we have to introduce higher order corrections to the the effective potential. Fixing 
our attention on the vertical motion only, and imposing ^^^ = 0, we obtain 

^y ~ ^yy"y + Ij^^yyy^y + '^^yyyy^y + ' ' ' (^L)) 

which clearly shows the coefficients connected to the anharmonicity. 
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To quantify anharmonic effects we use the ratio between the first nonlinear and 
linear term in the potential 

where the presence a, with a = A/L, reminds us that the importance of 
nonlinearities depends on the amplitude of the motion A. Using Equation (fTTni and 
changing the independent variables to x,Gy we can rewrite this quantity as 



L 

Pi = cr— 
9 
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(22) 



and its behavior for small oscillation frequency can be discussed looking at 
Figure (jH]). Near (a;*,G*) both terms go to zero, so the working point which minimizes 
the nonlinear effects is near the minimum of the curve. To first order, we must set 
Kyyyydy ^ Kyyy^S^ = 0, whcrc ^j arc the displacements from the critical point and K*, 
the expansion coefficients of the effective potential around it. Near a fold, however, only 
the first term goes to zero, and even worse, the second term increases without limits. So 
we expect huge nonlinear effects near a folding point, and a configuration near {x*,G*) 
should clearly be preferred. 

For finite values of frequency the more linear configuration does not coincide with 
the configuration that minimizes the frequency variation with the load. The reason is 
that a load variation changes both the working point and the inertia seen by the blade, 
and only the first effect is connected to nonlinearities. 

We conclude this section with an estimate of the importance of higher order 
nonlinearities. The ratio between the second order nonlinear term in the potential 
and the linear term gives 



^2 ^yyyy ^ 2 ^yyyy ( 9 



P2 = a^_^^a^^p _i^ (23) 

'yy ^y \^^y / 

which is dangerous because the numerator no longer goes to zero near the singularity. 

This result expresses the fact that when the vertical stiffness of the system goes to zero 

some kind of nonlinearity necessarily becomes dominant, and only the cubic term in the 

effective potential, proportional to its asymmetry, can be avoided. 

This can be seen as a lower limit on the vertical frequency obtainable, which we 

expect to be quite small in a concrete situation, because it is proportional to a^: 



/>" 



a K* 

y yyyy /o^^) 



27r\ L Gy 



As an indication, for a blade with L = 30 cm, a coefficient K*yyy ~ 10 estimated from 
a typical working point and an indicative amplitude for the motion A ~ 10^^ m ,we 
get / > 10~^ Hz, which seems quite reasonable. The coefficient Kyyyy is of course not 
the same for all the critical points, so it seems reasonable to find the best configuration 
from this point of view, if needed. 
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Figure 7. Sketch of geometric blade variations induced by a temperature change. The 
angle of the blade's tip is fixed due to the constraint of the the other blades in the 
MGAS. Li and L2 are the distances between each clamp and a fixed vertical reference 
axis placed in the center of the filter. Those distances change with the temperature 
according to the thermal expansion coefficients ai, a2 of the clamp materials. X is 
the horizontal length of the blade. 



2.1.5. Thermal Stability The temperature dependence of the blade equilibrium 
position depends on the mechanical and thermodynamical parameters of the material, 
i.e. the thermal expansion coefficient and the variation of Young's modulus with the 
temperature. A sketch showing how temperature affects the blade's geometry is shown 
in Figure 13 

Choosing the geometrical and mechanical parameters L\ ,L2 ,Q;i,a2 as shown in 
Figure m we have that x = X/L changes linearly with the temperature according to the 
following law 

x = xo{l + a^AT), (25) 

where a^, is an effective thermal expansion coefficient defined by the following formula 



a. 



^? 



LO 



— a, 



(26) 



with a the thermal expansion coefficient of the blade's material. For example, if we 
choose a = ai = a2 (supports and blade made out of the same material) we obtain 
C(x = 0. 

To first order, the material dependence on temperature enters through Young's 
modulus. In fact, we can write 

E = Eo{l + aEAT), (27) 

where a^ is a constant coefficient. Using the definition ^ we obtain the effective 
thermal expansion coefficient for Gf. 



ac = -{2a + oe). 



(28) 
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The Typical value for stainless steel expansion coefficient at ambient temperature is 
a ~ 1 X 10~^K~^. Common values for aE are negative and of the order 10~''-ft'~^[articolo 
NIM]. Accurate values strongly depend on temperature and on the alloy. This means 
that ac is positive and dominated by the Young's modulus variation coefficient. 

It is possible to tune ue both in amplitude and sign in a quite large range (say, 
aE G [—10~^K~^,10~^K~^]) by changing alloy composition and/or using thermal 
treatments [SMC 2003]. Such special alloys can be used indeed, to minimize temperature 
dependence of the working point. 

An interesting question is whether it is possible to find configurations which are 
insensitive to the temperature changes for a given set of coefficients a,. As our main 
issue is to obtain a system with a small vertical frequency we will limit our study to an 
understanding of the thermal behavior of monolithic blade near a good working point. 

We start with the dependence of the vertical position Y of the blade's tip as a 
function of the temperature. The first order variation of the vertical position due to the 
temperature is controlled by the quantity 

dY ^ , ^ dy , , 

which as we can see is the sum of two contributions. The first term accounts for the 
simple thermal expansion of the blade length. The second term is the "anomalous" 
contribution connected to the fact that the system is displaced from the working 
point by the temperature variation. Choosing as independent variables x and Gy 
(see [Appendix AD we can write 

dY f dy\ ( dy \ 

— = oY + Lx\ — \ax + LGy I -— - \aG = L{ya + '^^a^ + '^0^0) (30) 

The coefficient 7g can be rewritten as 

and it is quite dangerous, because it diverges near the critical points. This is very 
reasonable physically, because when the vertical stiffness of the blade becomes very small 
we can expect a strong sensitivity to the variation of blade parameters. One possible 
way to compensate for this effect is to properly tune the parameters in equation (|30| . 
We can adjust the term proportional to a by changing y, which is connected to the 
boundary conditions for the blade. As we will see in Subsection 12.1.61 there is a wide 
range of values for the y parameter, positive and negative, compatible with a critical 
point. 

The term proportional to a^, as discussed previously, can be tuned with an 
appropriate filter design. We can estimate its importance by looking at the left diagram 
in FigurelHl There the relevant quantity dy/dx = x'^j^ is plotted versus Gy for different 
X values in the stable region near the critical point {x*,G*) (compare with Figure IHj). 
It is evident that the rapid variation of this parameter around the critical point is 
particularly relevant when we are very near to it (x ~ 0.903, bold curve). This means 
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Figure 8. CoefEcients related to thermal stability. On the left side the quantity 
{dy/dx)Gy'is plotted as a function of the load parameter Gy, for different values of x 
choosen near the critical point. On the right side we plot uj~^{dujy/dx)Gy , as a function 
of the same parameters of the left side. See Section [2.1.51 for explanations. 



that it is surely possible to tune this parameter, but also that the tuning could be very 
delicate. 

Another important parameter, v^^hich should be made as much stable as possible, is 
the vertical oscillation frequency. We can define an eff'ective "expansion" coefficient for 
it in the following way: 



a. 



1 dujy 
ojy dT 






UJn 



' dujy 

dx 



GyOG ( dUy 



u„ 



9G. 



a 



(32) 



There are two quantities which depend on the working point in this expression. The 
first one is proportional to a^, which depends on the filter design, and to the variation 
of the logarithm of the vertical frequency with x, which we can rewrite as 



1 dUy 

ujy dx 



l_Gy_ 



dGydx ^ 



(33) 



This contribution depends fully on the nonlinearity of the system. The quantity (|33p 
is plotted on the right side of Figure |H1 as a function of Gy, for several x values chosen 
near the critical point. 
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The second working-point-dependent quantity is proportional to the variation of 
the frequency with the load. We met this quantity when we discussed nonlinear effects, 
and we saw that it can be set to zero by fixing the working point near one of the 
minima of Figure IHl Here the derivative is weighted with uj~^, so the situation is a bit 
different. As the minima of Uy are also minima for \ogUy it is always possible to set 
the contribution (|33| to zero. However, as we approach the critical point its variation 
will become very large. As for the vertical variation of the blade's tip, this means that 
the tuning of a^; will be possible but delicate, expecially very near to the critical point. 
The method used to evaluate all of the quantities that appears in these expressions is 
explained in [Appendix A 



2.1.6. Universality of the Results An important point is the robustness of the critical 
behavior of a monolithic blades. We want to understand what is the range for the 
"general design" parameters (tip angle, base angle and shape of the blade) which are 
compatible with the existence of a critical point. In order to do that we need a precise 
characterization of the "interesting" working point, that could be used to explore the 
parameter space. 

This characterization is provided by Figure where it is apparent that when 
{x,y) = {x*,y*) we have both 
/ flC \ 

{^\ y'') = K = (34) 



and 



, dy 



dy 



y_ 



x\yn=K,y = 0. (35) 



How can these conditions be extracted from Figure |31 on page [3? The curves there 
were obtained by continuing the solution in the x parameter, at a fixed Gy and with 
free y. We recognize that Kyy = when the slope of the depicted curves goes to 
infinity. This is a folding singularity that can be identified and continued adding a new 
free parameter, for example Gy, obtaining a locus of points where the Equation ^^ 
is satisfied. If, following this locus, we find another folding singularity, clearly there 
also Equation |r?K| becomes true. In this way the general strategy that can be used 
to explore the subspace of interesting working point is clear in principle: add another 
free parameter and continue the new fold. Our implementation of the method was a bit 
different because with AUTO2000 it is not possible to continue a fold of a locus of folds. 
It is based on the use of the auxiliary functions introduced in Section p\.ppendix A[ but 
in any case give us equivalent results. For each working point we know the value of all 
the parameters, and of course the solution for the equilibrium equations. Monitoring 
the auxiliary functions we can also determine to a given order the effective potential's 
coefficients K..., to get more quantitative information. 

Our preliminary tests show that there is a wide range of parameters where the 
MGAS concept is effective. This can be seen in Figure IHl which can be used to 
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understand the dependence of the working point on the blade base and tip angles. 
For a given pair of angles 9q^ 9i a working point exists only inside the admissible region 
delimited by the upper and lower bold boundaries. Inside the admissible region we can 
read the values of Gy and 1 — x by following the continuous and dotted contour lines. 
It is apparent that a very large range of possibilities exists, though probably very low 
values of Gy will not really be useful. 

Coming back to FigurelHon pageOwe note that the point {x*,y*) where the four 
bold branches with Gy ~ —1.6895 connect together is special. It is clearly a saddle 
point for Gy, and this tells us that the gradient of Gy is zero, which means 

K = (36) 

and also 

k;^ = K^ = (37) 

If the tangent to the curve at that point were vertical we could classify {x*,y*) as 
a so called "pitchfork" bifurcation point[Iooss and Joseph 1990, Arnold 1996], which is 
depicted schematically, with bold lines, in the upper right box. In the same box the 
dashed curve represents the locus of points where Kyy = 0, while the locus of points 
where Kyyy = is a horizontal continuous line which overlaps with the bold one. This 
means that the "interesting" working point is coincident with the "pitchfork" bifurcation 
point. 

However, a careful inspection reveals that what we really have a "two-sided 
bifurcation with a turning point". This one is depicted schematically in the lower right 
box of Figure HI emphasizing its different character. From this Figure we see that in 
this case the intersection between the dashed zero Kyy = curve and the continuous 
Kyyy = ouc docs uot correspond anymore to the singularity. More precisely we can 
expand Gy around the crossing point for the bold lines, obtaining up to the first order 
in 6j; 

^y ~ ^y = Q^yyyy"y + 2 [ vyy^ "" "^ ^yyy) "y + ^yyx^xOy (38) 

where in the pitchfork case (but not in ours) Ky^y = 0. Setting Gy — G* = we 
obtain the equation for the bold lines in the diagrams, while taking the first and the 
second derivative of Eq. (jHI^ with respect to y we find the equation for the dashed line 

yy 'n yyyy y """ \ yyy^^ """ yyy) y """ yyx"x U ['^^j 

and for the continuous one 

^yyy ~ ^yyyy^y "•" ^yyyx^x + ^yyy l^Uj 

which of course are valid in a small enough neighborhood around the crossing point. The 
main conclusion is that the real singularity associated with an "interesting" behaviour 
is a fold of the curve y{x, Uy = 0). As long as Kyyy is small (as it is in our example) we 
expect that this singularity will be near the crossing point. 
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Figure 9. Some possible working points for the system. The bottom angle is on 
the vertical axis, the top angle on the horizontal one. The set of working points that 
we found with the help of the continuation procedure are inside the upper and lower 
boundaries which correspond to the bold lines. For a given admissible pair of boundary 
angles 6*0 , Oi we can read the value of the vertical load Gy (on the continuous contour 
lines) and of the 1 — x parameter (the dotted contour levels). As an example, it is 
possible to obtain a good working point for a blade with 9o ~ 60 and 6*1 = by setting 
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A pitchfork singularity is the one which is usually used to describe Euler buckling 
of a beam under compression. In that case the buckling is associated clearly with 
the breaking of a symmetry which is not present in our case. But to completely 
understand in what sense the "interesting" working point we described can be associated 
to a "buckling" phenomena it would be better to have a connection with some kind of 
geometrical peculiarity of the blade's profile associated with the two different "buckled" 
configurations. Our conjecture, supported by some numerical calculations, is that in 
the bistable regime the two configurations differ in the presence of a point where the 
curvature changes its sign. 

2.1.7. Uniform Curvature At the working point the blade is bent to a profile 
determined by Equations (P) and (jZj). We are interested in configurations with optimally 
distributed stress. There are two main reasons for this. First, a blade with a uniform 
stress has a higher safety limit. Second, in the optimized case for a given load the 
maximal stress is reduced, which means that the intensity of creep and creak noise is 
also reduced (special alloys such as properly treated Maraging do not show an increase 
of creeping if the stress is increased). 

The dominant part of the stress is connected to the curvature radius R of the blade, 
so for an optimized configuration this quantity should be independent of ^. For the "old" 
GAS filter this is not possible[Bertolini et al 1999] because the boundary conditions 
imply zero curvature on the blade's tip. In the monolithic case we can substitute in 
the equilibrium equations (EI),© the uniform curvature solution 6 = Oo + ^{6i — do), 
obtaining an equation for 7(^) which can be explicitly solved: 

i9^-9oy-Gy[sm9{^)-smeo\-G^[cose{0-cos9oy ^ ' 

where 7(^) must be always positive in order to describe a physically acceptable 
solution. 

It is evident that the shape depends explicitly on the Gi parameters. This means 
that the blade profile should be tailored to the chosen working point, and it is probably 
very difficult to do this in practice. The solution we found tells us only that the uniform 
curvature solution is an equilibrium position of the system. To be practically useful this 
configuration must be stable, and with a low enough vertical resonance frequency. 

Care must be taken in writing the equations for the continuation procedure. The 
dependence of the blade shape from the parameters Gi should be considered in the 
equations described in [Appendix A| only to lowest order. In fact, higher order auxiliary 
functions ^(fc,/) describe physically the behavior around the working point of a well 
defined blade with fixed shape. 

We performed extensive numerical explorations of the constant curvature solution 
space, trying to find a critical point, but without success. It seems that the condition of 
constant curvature is not compatible with the low vertical stiffness property we want. 

After the discussion at the end of Subsection l2.1.6l this conclusion is very reasonable. 
We associated the critical point responsible for the vanishing stiffness with a transition 



Monolithic Geometric Anti-Spring Blades 22 

between two states which differ for the number of points where the curvature goes to zero. 
This means that near the critical point a constant curvature solution should approach 
a "flat" configuration, but this is possible only when 60 = di and the load is parallel to 
the blade's direction. This case can be studied analytically, and is incompatible with a 
vanishing stiffness. 

We take the negative result of our numerical experiments as a confirmation of our 
conjecture about the "buckling-type" nature of the critical point. 

3. Conclusion and Perspectives 

The main purpose of this article was a study of the low frequency behavior of a geometric 
antispring. We identified a singularity in the space of tunable parameters (load and 
geometrical constraints) which is associated with a vanishing vertical stiffness. 

Nonlinearities do not seem to be big problems for the kind of application we have in 
mind, where the ratio between the expected motion's amplitude and the blade's length 
L is very small. Cubic nonlinearities can be eliminated by tuning the working point. 

About thermal stability, our conclusion is that it is possible in principle to reduce 
the sensitivity of the system to temperature variations. The tuning can be done by 
carefully choosing the materials, the geometry and the working point. However, the 
quantities involved typically change very rapidly near the critical point, and this means 
that at the end the only possibility will be a case by case trial and error adjustment 
procedure, probably delicate and difficult. 

We have not yet performed an extensive exploration of the space of possible 
configurations. We suspect that this space is very large indeed, and the possibility 
of some optimization of the blade's behavior should not be neglected, if needed. Our 
next step will be the study of the behavior of the coefficients K*__ which describe the 
system around the critical point in this space. We plan also to study the system in the 
high frequency regime, including mass effect and internal modes. 

The numerical model described in the article has been implemented in a set of codes 
based on the continuation algorithms provided by the AUTO2000 library. This code 
will be used in the future to study and optimize special configurations, and is available 
on request. 
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Appendix A. First and second order response around the working point 

In principle information about the coefficients of the effective potential can be recovered 
by numerical differentiation of the series Gi{x, y) obtained with the continuation method. 
However, this method is very inaccurate especially if higher order coefficients are 
needed. For this reason we used a different approach to perform this task, which is 
basically a perturbative expansion of the external parameters around a given solution 
of Equation (jH]) and ^. 

There are many ways to choose independent parameters, and here we describe a 
possibility which is simple from a computational point of view. 

The more direct approach is to choose as independent parameters the external 
loads Gj, and to introduce a small perturbation e^ around a reference point by writing 
Gi = Gi + ei. For a given quantity C we introduce the notation Ci^m to identify the 
coefficients of its expansion in powers of e^. For example, we formally write the solutions 
of the perturbed static problem as 

0{i) = Y.^le%m- (A.l) 

k,l 

By direct substitution we can easily obtain up to a given order the equations for O^^i 
and pk^i. The equations for 6^0,0 and po,o are of course exactly the © and (0), with the 
same boundary conditions ^, and load parameters Gi = G\ . In the general case we 
can write 

^ = - 5R { [{iGf + Gf) r/., + ^r^,_i, + r^,,_i] e^^°'«} (A.2) 

^ = 7(0pm (A.3) 

where r]k,i are the expansion coefficient for e~*^'''°e*^ which are 

when /c, / > and zero otherwise. Note that the equation for pk,i depends linearly on 
Ok^i but nonlinearly on 6^1,1' with k' < k oi V < I. For example, at the first order we 
obtain 

^,0 = 'idifi, (A.5) 

'70,1 = ^^0,1, (A. 6) 
and at the second order 

r]2fi = id2,o - iflfl, (A.7) 

rii,i = ^^1,1 - ^0,1^1,0, (A.8) 

r?o.2=A2- 2^0,1- (A.9) 
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We also need boundary conditions for the newly generated differential equations. These 
can be simply obtained by imposing that the angles at the top and at the bottom do 
not change during the perturbation, obtaining simply 

0k,iW = 0kA^) = (A.IO) 

for fc > or / > 0. The coordinates of the blade's tip can also be expressed as a function 
of the 9k,i variables. By direct substitution in Equations (fTn|l and ^ we get 

x + ty = J2^'i f(^''°''Vk,id^ (A.ll) 

k,l •' 

Using these relations we can obtain all the physical predictions we need. For example, 

the effective stiffness matrix of the blade's tip (the Kij piece in the effective potential) is 

proportional to the coefficient which connects the first order variation of the coordinate 

with the load variation. Explicitly writing first order terms we get: 

ox f 

{K-^)xx = ^Q- = ^1,0 = - j ^1,0 sin ^0,0 d^ (A. 12) 

{K-%y = -^ = yo,i = J9o,icos9o,od^ (A.13) 

dx f 

{K-^)xy = -Q^ = 3^0,1 = - J Oo,i sin 00,0 d^ (A. 14) 

{K-')y. = -^ = yifl = / 01,0 cos 00,0 t/e (A.15) 

The symmetry in the indices is not apparent, but it can be demostrated, for example, 
assuming that the equations for 9k,i up to a given order can be obtained by minimizing 
a potential energy. 

An alternative useful expansion can be obtained by perturbing Gy = G^^^ + e^ with 
fixed X. Now the expansion of a generic quantity is 

C = Y.^[Q (A. 16) 

I 

and Gx is no more an external parameter, but a Lagrange multiplier which must 

be determined order by order to satisfy the constraint (fTn|) . 

Appendix B. Exact solution for constant width blade 

If the blade width is constant we can rewrite the potential in Eq. Q as 

where G = JGI + Gy, f3 = 9 — tp , cosz/' = ~Gx/G and sin-?/' = —Gy/G. This is 
formally the Lagrangian of a pendulum which evolves in the "time" ^, whose solution is 
exactly known in term of elliptic integrals. We can write the solution as 

sin^ = sn(^e + 0,A:) (B.2) 
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where sn(a:, k) is a Jacobi elliptic function and k,(j) two integration constants that can 
be determined by imposing the boundary conditions (A; > 0). It is easy to eliminate 0, 
and we get the condition 
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(B.3) 



where m is an arbitrary integer, K{k) the complete elliptic integral of the first kind and 
(Tj = ±1. We are interested in particular in the coordinates of the blade's tip, which can 
be written as 
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(B.4) 



where /(0, k) is the elliptic integral function of the second kind, and R{ip) the matrix 
which rotates a bidimensional vector by an angle ip. We also have to take into account 
that during a variation of the parameters the boundary conditions stay fixed. This 
gives us a relation, which can be obtained by imposing dV{\) = 0. Another constraint 
is that the coordinate x of the blade's tip does not change during the variation, which 
is equivalent to imposing dx{l) = 0. The last step is to evaluate the differential dy{l): 
in this way we get three linear relations between dy{l), dG^, dGy and dk that we can 
use to express 

dy{l) = xdGy (B.5) 

where x is a function of the working point, and x ~^ oo when the vertical stiffness of 
our costrained blade goes to zero. After some algebraic multiplication it is easy to see 
that we can express x as 
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(B.6) 



where all the quantities are evaluated at ^ = 1. In order for x to be divergent there 
are two possibilities: some of the partial derivatives inside the determinants could be 
divergent, or the 2x2 determinant could go to zero. Moreover, in the first case, as the 
determinant depends linearly on its elements, the only possibilities are a divergence of 
partial derivatives which appears in the 3x3 but not in the 2x2 determinant. 

We used the analytical solution for cross-checking purposes and as a convenient 
starting point for the continuation procedure. By introducing an interpolation 
parameter it is easy to continue the blade profile from the constant width case to the 
specific case we are interested in. 
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